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Abstract 

We present decoupled, separable forms of the linearized Einstein equations sourced by a 
string trailing behind an external quark moving through a thermal state of A/" = 4 super- 
Yang-Mills theory. We solve these equations in the approximations of large and small wave- 
numbers. 
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1 Introduction 



In [T], using the AdS/CFT correspondence [21 [3l H] (for a review see [5]), the expectation 
value (Tmn) of the stress tensor of TV = 4 super- Yang-Mills theory at finite temperature was 
calculated in the presence of an external quark moving a constant speed v. The calculation 
is based on the trailing string configuration of [HI E] (see also [S] for a closely related develop- 
ment, [9] for earlier work in a somewhat similar vein, and, for example, [10] for a brief survey 
of other recent literature on related topics). Conceptually, all that one requires in order to 
calculate (T^n) is a solution of the linearized Einstein equations in the AdS's-Schwarzschild 
background, sourced by the trailing string, with infalling boundary conditions imposed at 
the horizon and additional boundary conditions, appropriate to the absence of deformations 
of the lagrangian of A/" = 4 gauge theory, imposed at the boundary of AdS's-Schwarzschild. 
In practice, the calculation is somewhat tedious because there are fifteen linearized Einstein 
equations. In [T] these equations were somewhat streamlined from their original form in 
"axial" gauge, where all metric perturbations with an index in the direction orthogonal to 
the boundary are set to zero. However, the most interesting set of equations in [1] , namely 
the one that determines (Tqo), still involves four coupled second order differential equations 
for four unknown functions, together with three first order constraints. The complicated 
form of these equations makes them hard to work with, challenging to check independently, 
and difficult to extend to more general settings — for example, backgrounds where conformal 
invariance is broken or the motion of the quark is not uniform. 

There exists in the GR literature a highly developed "master field" formalism, dating 
back to such works as [TTl [T2l [T3] . whose aim is to obtain fully decoupled and separable 
forms of the linearized Einstein equations, as well equations for other fields with spin, in 
curved backgrounds with some symmetry. As a toy version of the calculation, consider a 
massive scalar, with action 



The first step is to consider a specific Fourier mode of and J in the x directions, with 





in the curved background 
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wave-number k. The equation of motion for this mode is 



^d^ + \dra^fdr - ^ 



m 



-J. 



(3) 



where k = \k\ and we have defined /(r) = h{r)a{ry. The corresponding master field and 
master equation are an equivalent way of writing ([3]): 



73/2, 



drfdr - Vrr. 



' Jry 



m 



k"^ 3 a', 3 a'2 ^ 3 a'' 



(4) 



2 a'' ' 4 a^ •' ' 2 a 
where primes denote d/dr and ^2 is the laplacian on the orbit spacetime 



dsl 



(5) 



The first aim of this paper is to reduce the linearized Einstein equations, sourced by the 
trailing string, to five master equations, each one similar to (jlj). The homogenous parts 
of these equations are well known, having been worked out in full in [T3] following earlier 
work including [151 [IS]- The source terms were partly worked out in [T7], whose methods we 
foUowE With an eye toward future applications to non-conformal backgrounds, we derive 
the master equations in the general warped background ([2]) dual to a finite-temperature 
field theory on R^'^, without making assumptions about the stress tensor that supports the 
five-dimensional curvature. 

With the master equations in hand, we specialize to the case of the trailing string in 
AciS's-Schwarzschild and reproduce some of the results of [1]. Also we consider an analytic 
approximation for the large k region, following methods developed in [11]. We explain in 
appendix |X] how the parameterization of metric perturbations used in this paper relates to 
the axial gauge parameterization used in [1]. 

The current work has some overlap with the independent study |20j . 



^Our methods are also reminiscent of those used in [18] to study sound waves in A/" = 2* gauge theory 
via a gauge-invariant formahsm. 
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2 Decoupling the linearized Einstein equations 



The first step in formulating tlie master equations is to expand the metric and stress tensor 
perturbations into Fourier harmonics involving scalar-, vector-, and tensor-valued functions 
on R^. As we describe in section 12.11 the scalar case involves nothing but plane waves, 
and the other two cases involve plane waves times appropriate polarization tensors. The 
coefficients of these harmonics are functions of t and r which in general have some tensor 
structure in the orbit space ([5]). We refer to terms in the Fourier expansions as scalar, vector, 
and tensor perturbations according to the type of harmonic involved rather than their tensor 
structure in the orbit space: indeed, one finds that the scalar perturbations are expressed in 
terms of a symmetric tensor in the orbit space; the vector perturbations can be expressed in 
terms of a vector in the orbit space; and the tensor perturbations can be expressed in terms 
of an orbit-space scalar. In the case of tensor perturbations, the orbit-space scalar is (up 
to an overall factor) the master field $7-, as described in section 12.31 In the case of vector 
perturbations, the orbit-space vector can be expressed in terms of an orbit-space scalar 
by use of one of the linearized Einstein equations, as described in section 12.41 Likewise, in 
the case of scalar perturbations, the orbit-space tensor can be expressed in terms of an orbit- 
space scalar $S) as described in section 12751 Each master field satisfies a second order master 
equation, which implies the linearized Einstein equations not already used in constructing the 
There is an intrinsic parity for tensor and vector perturbations which affects 
the polarization tensors but not the form of the master equation. Thus the five master fields 
are $5^^™, <l>5"'^'^,$y*'", and $5. The corresponding master equations (without reference 

to parity) are given in flM|l . fj^Tl) . and fl5^ . 

2.1 Scalar, vector, and tensor harmonics 

The defining equations for scalar, vector, and tensor harmonics on R'^ are 

(did' + fc')§ = (6) 
{did' + k^)Yj = d^Yj = (7) 
{d,d' + e)T,h = 9%-,, = = T^, (8) 

where indices are raised and lowered using the standard metric on R^. 

^We do not know of a general theorem that master fields satisfying second-order master equations can 
always be found. 
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The scalar harmonics can be chosen to be pure exponentials: 



§(fc,x) =e'*^-^. (9) 
They satisfy the normalization condition 



{S{k),S{k')) = {27ry6\k-k') (§i,§2)=/ d-^xSl{x)§2{x). (10) 

Vector harmonics may be assumed to be proportional to e'^'''^, but for each choice of 
wave-number k there are two distinct solutions, which are conventionally specified by their 
transformation property under the parity transformation x —x, k ^ —k, which preserves 
both the defining equation ([7]) and the exponential part of the harmonic. Denoting k = \k\ 
and k± = \/k\^ one may express the solutions as 

yf\k,x) = ^^e'^'' (q -ks k. 
These vector harmonics flTT]) satisfy the normalization conditions 

^ycvcn^^^^ycvcn^^/^)^ = {27rf6^{k-k') = ( V"'^*^ ( A?) , V"'^^ (fc') ) 



(12) 



under the inner product 



(Vi,V2) = / d'xY*{x)Yi{x). (13) 

JR3 



It is convenient for us to make the x^ direction privileged, because our eventual aim is to 
describe the emission from an external quark moving in the x^ direction through a thermal 
plasma of A/" = 4 gauge theory. A general vector field Xi{x) on admits the Fourier 
expansion 



X^(x) 



d^k 



Xr^k)Yr''{k, x) + Xr°(^)Vf "(Ai, x) + Xs{k)Siik, x) , (14) 



for some set of Fourier coefficients Xy^^{k), Xy'^'^{k), and Xs{k), where 

Uk,x) = -^dS{k,x). (15) 



5 



The expansion (fl^ is orthonormal with respect to the inner product (fT3!) . 

Tensor harmonics may also be assumed to be proportional to e^^'^, and there are again 
two distinct polarization tensors with definite parity: 



-kik2 



-kik2 -kihi 

U21.2 , 1,21,2 nl.2 , 1,2 

kokr 



-kiks kik'i 



kl 

~~k^ 



13 ^ 2 
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Trf^k.x) 



-/cs k2 \ 

2kik2k3, ki{kl — k^) 
P P 

kiikl — kl) 2kik2k3 



k' 



kl ) 



With the inner product 



2 Jr3 



we have the normalization condition 

A general symmetric tensor field Xij(x) on admits the Fourier decomposition 

d^k 



Xij{x) 



(27r)3 _ 



Xl{k)6ijS{k,x) +X^{k)Sij{k,x) 



where 



+ X^^™(A;)V?™'^(A;, x) + X°^'^(A;)V?f °(/e, x 
+ Xr"(fc)T^p(fc, x) + X°''''{k)T°f''{k, x) 



Y^pik, x) = —d^^Yf^ (fc, x) = -^5(,V°f 



k 

i{k,x) = -^didjS+ ^6, 



(16) 



(17) 



(19) 



(20) 



(21) 



and we use the notation {ij) = ^{ij+ji). The decomposition fl20l) is orthogonal with respect 
to the inner product ([T9l) . but not orthonormal because of some /c-independent factors arising 
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from the inner products of the derived harmonics in (12T!) . 



2.2 Fourier decomposition of the perturbations 

The background ([2]) is, by assumption, a solution of the five-dimensional Einstein equations. 



The explicit cosmological term is redundant because it could have been soaked into T^^] 
however, retaining it explicitly makes it easier to apply the formalism to an anti-de Sitter 
space example. T^i, includes any contributions from bulk scalar fields or other matter. Let 
the background metric and stress tensor be denoted G^^} and TjiJ . Consider a perturbation 
of both the metric and the stress tensor: 



where A is a formal expansion parameter which we eventually set to 1. 

Let indices h and c run over t and r, while indices i and j run over the three x directions. 
The metric perturbations can be decomposed into three functions h^c which are scalars on 
R^, two vector- valued functions h^i, and one tensor- valued function hij. These in turn may 
be expanded in Fourier series as indicated in f[T^ and fl20p : 




(22) 



G 



Gfl + \h 




(23) 




(24) 




(25) 



+ A^''=™'^(^, t, r) Vr°(^, x) + /,^'°<^'^(^, t, r) Nf'ik, x) 




- Hl{k, t, r)Sij S{k, x) + H^{k, t, r) Sij{k, x) 
+ i7^'^™'^(fc, t, r) '^(fc, x) + H^'°''''{k, t, r) V°f (^, x) (26) 
+ H^'^-^^k, t, r) Tip{k, x) + H^'^^^k, t, r) (^, x)] , 



where the factors of a(r) are chosen for convenience. Likewise one may expand 




(27) 
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Tbi{t,r,x) = a(r) 



(2^ 



T^{k,t,r)Si{k,x) 



+ r,^'^^°°(^, t, r)Vr-(^, x) + T^^'^^'ik, t, r)Yf''ik, x) 
d^k 



Kodd/ 



Tij{t,r,x) = 2a(r)' 



{k, t, r)6ij x) + (A;, t, r) Sij{k, x) 



It is also useful to express a general vector as 



(27r)^ 



;V^{k,t,r)S{k,x) 



Vi = a(r 



d^k 



y.odd/ 



rodd/ 



(27r)3 
+ Vy{k,t,r)Eii{k,x 

The diffeomorphism symmetry of the Einstein equations (122|) may be expressed as 



(28) 



(29) 



(30) 



(31) 



where v'^ is an arbitrary vector, and for convenience we have quantified the smallness of the 
coordinate deformation in terms of the same formal expansion parameter A that we used in 
the metric expansion (!23l) . Thus, at linearized level, (!3Tll becomes 



(32) 



For the purpose of formulating master equations for metric perturbations, we do not need 
to know details about how the zeroth order stress tensor arises from matter fields. Instead, 
we may extract Tj^J in terms of the unperturbed metric using the zeroth order Einstein 
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equations: 



7^(0) ^ p(o) _ IdWa^Co) , kMo) = HiflP--fT^°) t(°) t(°U 



= -A/-3V-^-//'-3-/ 
2 a a 

A a'2 3 a' /' 
Ti°) = Aa^ + /a'2 + 2aa7' + 2aa7 + \^ f . 
2.3 Master equation for the tensor modes 

Decoupled equations of motion for the tensor modes may be found directly by plugging the 
expansions fl26l) and fl29|l into the linearized Einstein equations. The result is the same for 
even and odd modes, so we will simply omit to specify parity in the following. The master 
equation is 

(□2_w)$t = -Jt (34) 

where 

$^ = c?l'^El Jt = 2c?l'^T^ 

, A;2 11a' , lla'2 \\o!' (^5) 
= 2A + - + --/' + -— / + --/ 
la 4 la 



,n^ A ^a!' 3 a'/ 
= - + 3— + ^— 



12 Q „i £1 (33) 



and, as before, ^2 = ^jd^ + drfdr is the laplacian on the two-dimensional orbit spacetime 



2.4 Master equation for the vector modes 

The results for vector modes are the same for even and odd parity, so we will not refer 
to parity explicitly in the rest of this section. As explained in [TTl [Tl], the first step is 
to consider quantities that are invariant under diffeomorphisms. To this end, consider the 
transformation properties of vector components of the metric and tensor perturbations, as 
can be derived by substituting the expansions (12^ - (130|) into (!32|) : 

fb fb + O'dbVv ^ - kvv 

t(o) ,^.(0) (36) 

V y I o V y rvJ- XX 

a a^ 
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One may form diffeomorphism-invariant combinations as follows: 



fV _ fV j_^f) ttV — \ zl^f) n 

Jb - Jb + i^'^b^T n - n + 



rp{0) 



If the are expressed in terms of the master field as follows 







A; 



V 



$y 



(37) 



(38) 



then the be and ij components of the linearized Einstein equations are satisfied automatically. 
The tj components lead to an equation for $y involving up to third order derivatives in 
r, while the rj components lead to a second order equation for $y. In deriving these 
equations, we found it convenient to pass to a gauge where Hj^ = 0, and to assume harmonic 
time dependence, e~*'^*, for f^, , f^, and $y (which are now indistinguishable from the 
corresponding unhatted quantities). A crucial point is that the third order equation for $y 
follows from the second order equation plus a relation that follows from the conservation of 
the stress tensor. More explicitly, from the 0{X) term in the equation 



(39) 



one may derive the relation in question: 



-D (a Tf, , 



+ 2T. 



(0) o ^(0) 

XX , l^r XX 



ka 



-T" + 



7,5/2 



V 



(40) 



where Di, is the covariant derivative with respect to the orbit spacetime ([5]). 

The master equation for the vector modes is just the second order equation following 
from the rj Einstein equations: 



[□2 -Vv) = -Jv 



(41) 



where 



Vv 
Jv 



4 a2 



2 a 



iV 



a 



t; + ^5.(/r^) 



(42) 



■^Two degrees of freedom, namely fY and can be reduced to one, namely $v, because of a first order 
constraint that would emerge if we substituted general into the linearized Einstein equations. This is 
explained, for example, in section IV.B.2 of [17] . 



10 



2.5 Master equation for the scalar modes 

For scalar modes, there is no longer a notion of parity. The diffeomorphism transformations 
of the metric perturbations are 

£ ^ + 2D^,v' f,' ^ - -vl + ad,v'y 



o Qj 

Because of the appearance of v'^ without derivatives in the transformation of and of Vy 
without derivatives in the transformation of H^, these quantities may be gauged away or 
used to construct diffeomorphism- invariant variants of /^^ and iff: 

M = hi + 2I^(.X.) Hi = Hi + \hI + ^fXr (44) 

3 a 



where 



^^ = l{f' + l9,Hl] (45) 



is invariant under diffeomorphisms with only Vy non-zero. Also, 

- + v'^^T^^ + 2Ti°)a.).^'^ 



and the gauge-invariant quantities are 



1 

2''' XX ^^T 



(47) 



In manipulating the linearized Einstein equations, we employed a gauge where = 
and Hi = 0. One may however pass immediately to diffeomorphism-invariant formulas by 
replacing fl^, Hi, , , and r^^ by their diffeomorphism-invariant relatives as defined in 
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and (H7|) . and we will quote all formulas in diffeomorphism-invariant form. One may 
eliminate fff by noting that the ij Einstein equations with i ^ j — that is, the off-diagonal 
equations in the lower-right 3x3 block — lead directly to 

frS _ ^ nbc fS _ ^" -.s ('48) 

The tr, rr, ti, and ri Einstein equations involve only first order derivatives of and using 
stress tensor conservation equations one may show that the other Einstein equations are 
implied by these four. Approximately following [13] we express the symmetric tensor /^^ as 

M = l(2X-Y) ff^^'^Z f^ = ^i-2Y + X). (49) 
Then the first order equations just mentioned reduce to 

X' X+ [2 ^]Y + -\ — - + — ^ ] Z = Tx 



Y' - —X + + —Z 




2/ 2f P ' (50) 
^^ + V + ^/ + ?-//'V+P + C-^/-^/-30^-3-//')F 



43 2a / \ 4a^3 a 



where 



2 . 2a2 . 2a2 / ^ /' A . 

-af^^ - —f^ + — 2adr - 3a- + 10a' 

'UJ rv K \ J I 



TX 

Ux) 

2a^ ^ c 8a^ , „ r> /\ s 
ty = -r'Tr + -JT + 2a 

2a2 5 Sa^ o 



(51) 



rc = --^r, + 2a/r,,-^r, + ^ ^12.. +3/ -—f-^-jr . 

Because the system of equations ( l50l) consists of three first order differential equations and 
one algebraic constraint, it can be reduced to a single second order differential equation 
satisfied by any desired linear combination of X, F, and Z. In order to put this equation in 
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a form that resembles the other master equations, we define the scalar master field as 



.a' 



<^s = A\^X + Y + 3-Z ) , (52) 
where 

-"- H^'HHt'Lfa'r^ //.P.2Aa^.3/a« + £ar. (53) 
The master equation then reads: 



{n2+Wsdr-Vs) ^S = -JS 



(54) 



with ^2 = jLu'^ + drfdr and 



Ws = ^ (4Aaa' + 9f'a'^ + Afa'a" + 5a fa" + Saa'f") . 



(55) 



In the general warped background ([2]), the expressions for Vs and Js are too long to be 
reproduced here. Starting in section [31 we specialize to the case where the unperturbed 
metric ([2]) is AdS's-Schwarzschild, namely a = r/L and h = l — rj^/r^, where A = — G/L^. In 
this case Ws = 0, so the master equation fl5^ takes the standard form (^2 — Vs)^s = ~Js 
where 



Vs 



Js 



108rf^(5r^ + 3r|^) 



/ 



a^ (lAr% 2r'^ + 4:rjj 



( 72rjj 3(r^ + hr%) H 



14^2 



f 



(56) 



(Note that H as quoted in ( l56il is the same one defined in (l53ll .) 



3 The trailing string in A(iS'5-Schwarzschild 

Now, as promised, we turn to an application of the master field formalism to computa- 
tions of (Tmn) as sourced by the trailing string configuration in ArfS's-Schwarzschild. In 
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section 13.11 we work out the stress tensor in five dimensions created by the traihng string, in 
the Fourier-expanded forms (127|) - (|2U]) . The resulting expressions may be inserted into the 
master equations flM|l . fHTl) . and fl5^ to obtain explicit differential equations in r for the ra- 
dial dependence of the master fields, which we give explicitly in fl68|) - fl72|) . In sections I3.3H3. 51 
we obtain asymptotic forms near the horizon and near the boundary of AdS's-Schwarzschild 
for the solutions of the master equations. As explained in section 13.21 the asymptotic forms 
near the boundary translate into the desired quantities (T^n) in the gauge theory. 

3.1 The stress tensor of the traiUng string 

In Poincare coordinates, the 74(iS'5-Schwarzschild metric can be realized as a particular case 
of ([2]) with 

a(r) = y h{r) = 1 - ^ • (57) 

Here, denotes the location of the black hole horizon, and the AdS radius L is related to 
the cosmological constant term via 

A = -^. (58) 

It is useful to note that horizon temperature is T = r/z/vrL^. The equations of motion follow 
from the action 



(59) 



where Qajs is the worldsheet metric. In static gauge, the worldsheet metric that minimizes 
the action is, to leading order, given by [6l [7] 



9ap = G^,d^X^dpX'' X^{t,r) = {t r X\t,r) o' 
X\t,r)=vt + ar) ar) = -p^Ug^-^^ + z\og^-±^] . ^^^^ 

We can use the above expressions to compute the source term T^^ that appears in the Einstein 
equations (1221) . Because AdS's-Schwarzschild is the solution to the Einstein equations with 
negative cosmological constant but vanishing stress-energy tensor of the matter fields, T^^ 
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comes only from the trailing string: that is, 



rpjjLU 



2'Ka' 



2 7-0 











V 





o\ 








L^h 


h 








[v^ - h) 


rfjV 














L2 




V 

h 




rfjV 
L2 





























\ 
















0/ 



1. The notation 6^(x'' 



(61) 
X') is 



where we have set the formal expansion parameter A 
short for 6 (x^ ~ vt — ^(r)) 6{x^)6{x^). The expression (l6Ti) (with indices lowered) can then 
be decomposed as in (I27|) - (l29|) . with the coefficients of the decomposition given by 



y,odd _ V,odd 
'a ' 



T 



T.odd 



r 



_V,cvcn 



= T 

r 



6L2 



T 



r 







3kl 



-T 



krh [rj,/L^ 



tr 



r 



h 




T 



V.oven^ 



V.cvcn 



Ak^L^ 



„T,even 




rt 



T 



h-v^hy/L^ hv^rjj/L^ 
hv'^r'jj/L'^ v"^ — h 



(62) 
(63) 

(64) 
(65) 



where the common prefactor T is 



27ra' v^T^ 



(66) 



The odd components are zero because of the choice of the even and odd vector and tensor 
harmonics that we made in section 12. 1[ As a result, the odd-parity master equations may 
be solved trivially by setting ^'^'^ = ^'^'^ = 0. We will henceforth be concerned only with 
the even parity cases. 

The time- dependence of the source is e"™*^!*, and we will make the steady-state assump- 
tion that the master fields $t, ^v, and $5 have a similar time dependence. It is convenient 
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to make the following definitions: 

-- 

$5 = 



2iKlvkiL^ 
QkIL 



-ivkit 



■ivk\t 



(67) 



The equations satisfied by ijji with I = T, V, and S can be written as 



1 



drtt^fdr + 



v'kf k' 4 - uj 



f 



L2 



Vi 



-Jl 



(68) 



where i^t = 2, i^y = 1, and z/5 = (compare with equation ([3])). The quantities V/ and J/ 
are given by 



Vt = Q 



9r%_ 



12r| 



2k^ 
1^2 



(69) 



k'^L'^r^ 



Jv 



k±e 



-iki 



kk\ 



vki 



5vkirjj 



Js 



(70) 
(71) 



ki {k\2 + _ 3^2^2) ( -3A;2l^1 + [k''L\2k^ - ^v'^kl) - 90r^) ^ 



SzveL'^rl {2k^ + t;2A;2(p _ 3A;2)) J_ + 3A;iLV|, (l8i;^fcf + k\2 - 5v' 
3v^ekl{l - 2v')) 1 + Qivr%L' {2k* + 9v^kl{e - 3kl)) \ 



+ 18kir% (fc^(2 + 3v^) - 9v^kl) 

(72) 

where H = k"^ + 6rjj/L'^r^ as in 0561) . Note that the Vj decay at the boundary as l/r2 or 
faster. 



3.2 The holographic stress tensor 

To understand the holographic relation of the master fields to the gauge theory stress tensor 
Tmn, a useful preliminary is to count the independent components of Tm„. Because the 
stress tensor is symmetric, we start with 10 components. Conservation eliminates four, and 
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confer mal symmetry eliminates one more, so there are five left. These five independent 
components are related holographically to the five master fields. One can be more explicit 
by expanding 

/d^k 
j^^qs{k,t)S{k,x) (73) 



(Tor) conserved / /cy ys 



(2vr) 

(-^ij) conserved 2 



^dtqsik, m + qr^'ik, t)Yr''{k, x) + qf\k, t)Yf\k, x) 



(74) 



d^k 



^^^6,,Sik, x) + -^{38^ + e)qs{k, tp,,ik, x) 



(27r) 

tmr^''^ ^) + Istgy'^ik, t)v°f (^, x) (rs) 



+ t)T?f °(fc, x) + q°T t)T'^j{k, X 

The first term on the right hand side of (I71|) must take the form indicated in order to 
have 5'"(Tmo) conserved = 0. The first term on the right hand side of (ITS]) is dictated by the 
conformal condition (T"^m)conserved = 0, and the next three terms are consequences of the 
conservation equations 9™ (Tmi) conserved = 0. 

The boundary asymptotics of each master field encodes the value of the corresponding q 
coefficient in fl75]) - fl75]) . as we explain more precisely in sections [3.3H3.5I We have already 
remarked that the odd perturbations vanish, and so in the following we omit parity labels. A 
subtlety in the expansions (173!) - (fr5l) is that they actually don't capture the full stress tensor, 
which is 

{^mn) (-^mn)bath ^~ (-^mn) conserved ~l~ (-^m7i)drag • ('''6) 

Here (Tmn)bath = ^^^^^^ diag{3, 1, 1, 1} is the stress-energy dual to the black hole in AdS^- 
Schwarzschild. The last term is not conserved and is due to an external force required to 
balance the drag force of [6l [7]. The presence of this term was noted in [1], and we rederive 
it in section 13. 6[ 

In order to compute (T^n) using AdS/CFT, the most straightforward approach is to 
express the perturbation fl23l) in axial gauge, hj-fj, = 0. As will be explained in detail below, 
the non-zero components of the metric perturbations have the following large r behaviorj^ 

hmn = RmnT^ + Pmn' + Q^nk + (^(^"') , (^7) 



"'The quantities Pmn, Q^mm ^'iid Rmn differ from the corresponding quantities in [T] by different normal- 
ization factors. 
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where Rmn, Pmn, and Q^J^ may depend on x"^ = (t, x) but not on r. Up to a subtlety regard- 
ing residual gauge freedom]^ the boundary condition Rmn = corresponds to requiring that 
the metric on the boundary should remain fixed, or in other words not adding components 
of Tmn to the gauge theory lagrangian. Another boundary condition comes from requiring 
purely infalling modes at the black hole horizon. 

Because h^n is a perturbation around Ad^s-Schwarzschild rather than pure AdS^, it 
relates only to the second and third terms in (176!) . More precisely, as argued in [1], one has 

(-^mn)tot = (-^mn ) conserved ~l~ (-^mn)drag 27-3 Qmn ' C^^) 

The notation {Tmn)tot is meant to indicate the total VEV arising from the trailing string: 
it excludes (Tm„)bath- Actually, it excludes one other thing: the Pmn/f term in (177|) signals 
the need for a counter-term subtraction to make (Tmn) well-defined. As discussed in [1], this 
subtraction is a delta-function supported at the location of the quark. 

The main purpose of sections I3.3H3.5I is to relate the boundary asymptotics of solutions 
to the radial master equations fl68l) to the quantities Q^^ entering fl77|) . In so doing we will 
discover a natural split of Q^n i^^o a non-conserved piece which can be computed analytically 
and a conserved piece which must in general be obtained by solving the master equations 
with appropriate boundary conditions. We also explain what these boundary conditions 
are, thereby completing, in the master field formalism, a specification of the boundary value 
problem that determines the holographic stress tensor. 



3.3 The asymptotics of the tensor master equation 

Near the boundary of Ad^s-Schwarzschild, the solution to the tensor master equation behaves 
as 

Mr) =Rr{l-^,-^,logry Pt^ + Q^^, + O(r-) , (79) 
where Pt and a are given by 

Pt = |^ a^L\e-vX) (80) 

^The subtlety is that the boundary condition Rmn = is actually a little more restrictive than necessary. 
It should be allowed to have Rmn oc gmn, corresponding to adding T™™ = to the lagrangian, or Rmn oc 
V{mkn)^ corresponding to adding a multiple of dmT"^^ to the lagrangian. Such deformations correspond to 
metric perturbations in AdS's-Schwarzschild which are pure gauge. 
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and and Rj- are the two integration constants. Since H^'*^^™ = ^hilL^ipTir)e the 
requirement Rmn = translates into Rt = 0. Assuming Rt = 0, we have 

which can be used to find the tensor components of the Q^^^ coefficients via (12^ - (126|) . 
Feeding the result through (I78p and comparing to the expansion (175]) . one finds that 

gT(fc,t) =QT*(^)e-™'^*- (82) 

(Recall that because odd parity perturbations vanish for the trailing string, we are always 
referring to even parity; however, for a more complicated gravitational source, fl82l) could be 
used equally in reference to odd parity modes.) The normalization factor between $t and 
ipT in fl67j) was chosen to make the relation fl82l) simple. 
Near the horizon, the leading behavior of ipx is 

^^(r) = UT{r- ^^)-^-^i^V4r^ + VT{r- r^)™'=i^V4rH ^ _ ^ (83) 
and the infalling boundary condition is = 0. 

3.4 The asymptotics of the vector master equation 

In the vector case, the master field has the following asymptotic behavior at large r: 

iJv{r) = Rv- (l - logr) + Py-^ + 4 + ^(^"') ' ^^^^ 

where 

p^^J_v_vJ^ ^^L\e-v^kl) (85) 

Kg K 

and Qy* together with Ry are integration constants. Relating this large r series expansion 
to (177|) is not hard: one can use the definition (167|) and plug (18^ into (!38|) to obtain a large 
r asymptotic expressions for the gauge-invariant quantities /{,, and then use the expression 
for fh in (1371) as well as the axial gauge condition = to find large r series expansions 
for both fY and . When the latter quantities are plugged into one can easily 

check that all non-zero components of Rmn turn out to be proportional to Ry, so setting 
Rmn = means setting Ry = 0. 

When Ry = 0, the large r behavior of the functions that enter in the metric perturbations 
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is given in axial gauge by 



(86) 



Feeding ( 186|) through (ITHl) and comparing to (1741) . one finds that 

qvik,t) = Q'^\k)e-'^'^^\ (87) 



provided one identifies 

<^«->--|f (88) 

At the horizon, the situation is similar to the one encountered in the tensor case. The 
leading behavior of ipv is 

^y(r) = Uv{r- r^)"-'=i^V4r-« + Vy{r- thT^^''"'''" + 0{r - m) , (89) 

and again the infalling boundary condition is VV = 0. 



3.5 The asymptotics of the scalar master equation 

Near the boundary, the asymptotic behavior of ips{r) obeying fl68l) is 

Mr) = + Rs^ logr + 1 + 0(r-3) , (90) 
with Ps explicitly given by 

' k^i2 + v^)-3v^kl 



71 

^5 

itot 



Ps = -, ' ' \ (91) 

Ki 2k^ 



and Q^* and Rs being regarded as the two integration constants of the second order equation 



The asymptotic behavior of ipsi^) given above can be related to the asymptotic behavior 
of the metric perturbations fl77|) through a straightforward but tedious calculation: one just 
has to use fl67|) to go back to the master field $5 and then trace back the steps outlined in 
section 12.51 in the special case where the background is AdSs-Schwarzschild. It is important 
to note that in doing so we only need to consider large r series solutions to the differential 
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equations we encounter — a much easier task than solving these equations for all r. We will 
now provide an outline of the computation. 

First, one can easily obtain expressions for /^^ and in terms of il)s{f) that are valid for 
all r from combining the definitions (H8|) and fHOj) with equations fl50l) and with the definition 
fl52|) of the master field. Next, one can find large r asymptotic expansions for /^^ and Hf^ 
by replacing ips{i^) with a large r series solution to the master equation. While the first few 
terms of this solution are given in (pOl) . an improved series expansion (up to order O 
is needed to arrive at the results listed below. The last step of the computation is to assume 
axial gauge, namely = f^^ = = 0, and to use the asymptotic expansions for /^^ and 
Hi to find series solutions for the quantities /J^, ff, H^,, and iff that obey the system 
of differential equations (jH]). The asymptotics for the non-zero metric perturbations hmn 
can then be found from fl2^ - fl26|) . and one can check they have the form (1771) . Moreover, 
all scalar components of Rmn end up being proportional to Rs, so again Rmn = means 
Rs = 0. 

Assuming Rs = 0, the expressions one ends up with for the axial gauge metric pertur- 
bations are 



9 r 2L V 4kiL^ J 



Hi . ,.-e-^- + ^ ^-gr + j e--- + (92) 

~ ^ r3 + 2L V 4P 

As in the tensor and vector cases, feeding (!92|) through (178|) and comparing to (173ll . one finds 
that 

qs{k,t)=Qf{k)e-'''^\ (93) 

provided one identifies 

(r„),„..|!^/|^jle-*'S(*,2) (94) 
(T,,),„. ^ 2|!^ / (gi^%S(.:f) + ^^S,(«.f)) . (95) 
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From (p2l) one may also find explicit expressions for the quantities Pmn appearing in (177|) . 
Using fll42p in appendix |X] one can check that the expressions so derived match with the 
corresponding Pmn coefficients calculated in [1]. 
At the horizon, the series expansion for ips is 

Mr) =Us{r- r^)-»^i^V4.« + _ ^^^-k^LV^r, ^ (j^^ _ _ (gg) 

The inf ailing boundary condition is = 0. 



3.6 The force on the moving quark 

In sections [3^ and [3751 we saw the non-conserved part of the stress tensor, ^mn) Ats.^-, emerg- 
ing from mismatches between the Fourier expansion (I73p - (175]) of a general conserved, trace- 
less stress tensor and the boundary asymptotics of the metric perturbation /imn, as computed 
by tracing backward through the master field definitions in sections 12.41 and 12.51 In this sec- 
tion we recheck an observation of |T], namely that the non-conservation of (Tm„)drag is the 
result of an external force on the quark that precisely cancels the drag force so as to keep 
the quark moving at constant velocity. Equivalently, one can regard this non-conservation 
of (Tm,n)drag as describing the rate at which energy and momentum are supplied to the bath 
by the external quark. 

The external force on the quark is 

r=\ ci\9„(T™)d,ag, (97) 
Jv 

where the integration volume V is any three-dimensional volume enclosing the quark. The 
quantity <9„^(T™") can be easily calculated from (IMl). and (1^: 

^5 



dml(r-^-^) = ^^b{%^ -vt)b{x')b{x'')[v 1 0). (98) 



Evaluating the integral and using = — f^, we obtain 

which is indeed minus the drag force P, E]. In obtaining the last equality we have used 
T = r/^/vrL^ as well as L"^ /a' = ^/g^^J^ . 
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4 Small k approximation 



If /c <^ T, the equations for the ipi simphfy and can be solved analytically. In fact, we can 
use the regular perturbation theory technique to find small k series solutions for each of these 
equations. The series solutions will in turn enable us to extract small k series expansions for 
Qt*' Qv^y which are needed to find the small k behavior of {Tmn)- 



4.1 Tensor set 



As remarked before, the tensor equation is the same as the dilaton equation in [21], or the 
A-equation in [1]. While the derivation of the small k solution of this equation can be also 
found in section 3.5 of [T], we will now outline the main steps of the calculation. 
We start with the equation for ipT in (l68il and develop a series expansion in k: 



iprir) = ipoir) + k'ipi{r) + k'^'ip2{r) + . . . . 



(100) 



The equations satisfied by each ipi can be found from equating like powers of k in fl68p . In 
doing so, we consider ki and k± to be of the same order of smallness as k, so for example a 
term containing kik^k"^ will be considered to be 0{k^). These equations can then be solved 
by requiring that the boundary and horizon asymptotics of the solutions be nothing but 
small k expansions of fl7^ and fl55]) with Rt = and = 0. For example, the equation for 
V'o is 

^0 = -^^, (101) 



hi\ k'^rh 



and its solution is 



Mr) = — ,, 2 3 2 tan vr + log 



log- 



H 



Expanding this solution in series at large r, we obtain 

4 v^kl 1 4, v'^klm 1 



0[r' 



(102) 



(103) 



from which we can identify the coefficient of 1/r with the leading contribution to 



tot 



23 



Performing a similar analysis for ipi and 1/J2, "we find 



/^tot 



4P 



8P 



+ 



t;^A;iL4(67r- 121og2) 
192A;2rH 



2 ,,2 



TT^ - 12(log2)^ 
67r - 12 log 2 



(104) 



4.2 Vector set 

The small k analysis in the vector goes mostly the same way as in the tensor case: one starts 
by expanding ipy in a low-k series 



il;y(r) = ^o{r) + k-^iir) + k'^'^2ir) + ... , (105) 

and then the equations for ipi follow from equating like powers of k in the expansion of the 
ipv equation in fl68l) . The difference from the tensor case is that the matching of the tpi to the 
asymptotics fl8^ and fl89|) with Ry = Vy = and appropriately expanded in k results in one 
undetermined integration constant at each order i. Fortunately, this integration constant 
can be determined simply by looking at the horizon asymptotics of For example, ipQ is 
given exactly by 

V'o W = yr, + f^o^ , (106) 

where Uq is an integration constant of order 0{k^). It is not hard to see that, to order 0{k^), 
this expression satisfies the required boundary conditions: ipQ looks like an infalling wave at 
r = th and there's no 1/r piece at large r regardless of the value of Uq. The only way to 
determine Uq is to go to order 0{k). 



We obtain the following small k expansion for Q^"*" 



^ 4 ruk^i^v'kl - k^) zk^L' [k^ - (81og2)^^fct + {2\og2yekl] , 
kI Avkkl kI Wkkf ^ ' 



(107) 



4.3 Scalar set 

Obtaining a small k solution for the scalar equation in fl68|l is a bit tricky, because the quantity 
ips doesn't have a small k expansion itself. Instead, we can find a small k approximation for 
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Mr) = ik' + i^,)Mr) . (108) 



the rescaled field 

6rj 

To do so, we can follow the same steps as the ones outlined in the tensor case in section |4?T] 
and replace ipT by ips- It is not hard to check that the zeroth order solution ipo is zero in 
this case, so the first non-trivial term should be first order in k: 

ipsir) = kipi{r) + k^ip2{r) + k^ipsir) + ... . (109) 

One has to be careful though when matching the solutions for ipi to the asymptotics (1901) 
and fl96l) : the rescaling factor in (llOSp reduces to k^ in the large r limit, so the order of the 
ips expansion changes by two when we look at the boundary asymptotics. With this minor 
technical difference, the matching proceeds in the same way as in the tensor case. 

By reading off the coefficient of the term in the large r expansion oiips^ we find Qg* 
to be 

^ L tvrUk' + ^kl) 3v'rHkn-^v'kf + k\2 + v')] 

4L^kiik^-3v^kl)^ 4 {k^-3v^kiy l^^uj 

This formula diverges at the Mach angle, namely where the ratio of the momenta is ki/k = 
l/v^/S, signaling a sonic boom in the thermal medium. While the divergence appears to get 
worse as one goes to higher orders in the expansion, the singular behavior disappears if the 
series expansion flllOp is resummed into a form that resembles a Lorentzian lineshape. For 
a more detailed discussion, see section 3.5 of [1]. 



5 Large k approximation 

In [in], a method of Green's functions was employed to determine the large k asymptotic 
form of the VEV of the operator dual to the dilaton, as explored numerically in |21j . 
This method can be understood as an expansion in powers of rn, which makes sense because 
th ccT and all /c-dependence of dimensionless quantities arises in series expansions in k/T. 
Suppose one starts with an ordinary differential equation of the form 

£^(r) = [dl + t(r)(9, + n(r)] ^/^(r) = -]{r) (111) 
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where there is an expansion 

C = Cq + thCi + r]jC2 + ■■■ 

iJ = tpo + TH^i + rjjiJ2 + . . . (112) 
3 = jo + rnji + rHj2 + ... ■ 

The radial master equations fl68l) fit this form after both sides are divided by / = (1 — 
rjj/r^)r'^/L'^. At least formally, one may solve (11 111) order-by- order in rn. At the lowest 
order, Cq is a differential operator in pure AdS^, which is the th ^ limit of AdS^- 
Schwarzschild. Construct two solutions of the homogeneous equation. 



where iJjh satisfies appropriate boundary conditions at the degenerate horizon r = 0, while 
ipYEY has appropriate falloff at the boundary of AdS^ to describe a VEV rather than a 
deformation. Then a Green's function inverse to Co, satisfying 



as well as the boundary conditions at the horizon and infinity satisfied by ipn and ipyEY, 
respectively, may be constructed as 



ColpH = = CotfjYEY , 



(113) 



£oGo(r;ro) = S{r - tq) 



(114) 



Go(r;ro) 



1 \i^H{ro)ijYEy{r) for r > tq 



(115) 



l^H, ^YEY]{ro) UvEv(ro)^/'H(r) for r < tq 



where the Wronskian is 



W[ijH, V'VEv](ro) = ^H(ro)^vEv(^o) - ^vev(^o)V^h('^o) • 



(116) 



The leading order solution i/jq may be determined as follows: 




(117) 



where 




(118) 
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In optimal circumstances, Bq{oo) = lim^^oo -Bo(^) is finite and Ao(r) fast enough as 
r — > so that the boundary asymptotics of ipoi^r) may be expressed as 

^o{r) = Bo{oo)^pYEvir) + (sub-leading) . (119) 

The form flll9p is convenient because, up to a normalizing factor, Bq{oo) is just the 0{r^) 
contribution to Q^^, Qy^, or Q^g^. Unfortunately, there is a complication in the cases we 
will study: BQ{r) does not have a finite limit. The solution is to isolate a simple analytic 
form ipQ^^lr), capturing that part of ipoir) whose boundary asymptotics are more singular 
than indicated in (11191) . and split both tpQ and jo so that 

^0 = ^^^"^ + ^^ Jo = Jo'^^ + Jo'' 

r ,/,CT _ _ -CT ,R _ _.-R ^ ' 

'-oV^o — Jo '-oVo — Jo • 

One may then follow through the Green's function analysis of the "renormalized" equation 
£oV^o' ~ ~Jo' indicated in flll7l) - flll8p and find that the "renormalized" analog of (11191) 
holds with finite B^{oo). The splitting (fT20|) is in the spirit of holographic renormalization 
as reviewed for example in [22], although it also follows an older and more elementary line 
of thinking where one looks for leading non-analytic contributions to Green's functions in 
momentum space: see for example [3]. 

Higher-order ijji may be determined in a similar fashion: 

i^iir) = -{Go* ji){r) j^=j^ + CiiljQ 

^2ir) = -{Go* j2){r) 32 = 32 + C2ipQ + Ciipi (121) 



A "renormalization" process similar to (11201) could be employed, if needed, at higher orders; 
but need for it would be unexpected given that the "counter-term" contributions usually 
correspond to contact terms in the gauge theory which are temperature-independent. 

A significant obstacle to the treatment just described is that the expansion (11121) is not 
uniform in r: near the horizon it breaks down. It is not clear whether the boundary condi- 
tions appropriate to the zero-temperature Green's function Go lead, through the expansion 
(11121) . to a solution if) satisfying standard boundary conditions at the horizon. This issue 
was confronted in [19] with a WKB treatment. We will take the less rigorous approach 
of computing the leading non-trivial /c-dependence of Q^*, Qy*, and Q^* and comparing 
to numerical evaluations from [T|. The eventual aim is to give an all-scales description of 
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(Tmn)-, both in k space and position space, using a combination of numerical methods for 
intermediate scales and analytical approximations for the IR and UV limits. 



5.1 Tensor set 



This is essentially the case considered in [19] , because the tensor master equation is the same 
as the dilaton equation [21]. The leading non-zero terms in the expansions (1112 1) are 



Jo 



J2 = —^Pt 



(122) 



where Pt is as defined in f lHOj) and 



(123) 



The zero temperature solutions and Wronskian are 



1 



1 



^VEv(r) = ^h{kL'/r) ^^(r) = -K^ikL'/r) 



iy(r) = iy[^H,V'VEv](r) = -^ 



(124) 



where Ip and denote Bessel functions. It's significant that fc^ > always because it 
results in a Green's function that vanishes faster than any power of r as r — 0, i.e. in the 
depths of AdS^^ Using ^/'^'^(r) = Pt/t^-, one may straightforwardly show that 



Mr) = Pt 



SnkL^ 
16r4 



0(r- 



ivki 



(125) 



A "counter-term" is not needed for the analysis of 1^2(1^), as expected on general grounds. 
Comparing (11251) with (!79l) leads immediately to 



Q 



tot 
T 



T 



SnkL"^ ivki o 4 n 

+ ^rjj + O(r^) 
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PL2 



(126) 



^Probably this means that putting the horizon back at finite corresponds to 0{e ^^'^/^h^ effects for 
some c- number 60. If so, it means that the series expansion (jll2|) for in powers of rn is only an asymptotic 
series. Because th/L^ = t^T, one sees that e"'^^ has essentially the form of a Maxwell-Boltzman factor. 
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5.2 Vector set 

The leading non-zero terms in the expansion flll2p are 



Jo = J 2 



(3F - v'kt)Pv 



where Py is as defined in ( 185|) . The zero temperature solutions and Wronskian are 

W{r) = W[ijH,i^YEv]{r) = -^. 
Using ipQ^{r) = Py/r'^ one obtains 



Mr) = Pv 



1 TTkL^ 



47-3 



Mr) = Pv 



3F - v'kj 
Sivkik'^L'^r'^ 



0{r- 



Comparing to one sees that 



QT = Pv 



TikL' 3k'-v'kl 2 



4 Sivkik'^L'^ 



r'n + 0{t],) 



5.3 Scalar set 

The leading non-zero terms in the expansion flll2p are 



£0 = dt + -dr -—k' + - 
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1 ^ {2k'' - v^kDL" 



3r^ 



2iL'v 3k^ - bv'kfil - v') + kfk'{l + 8v') 1 tvk^ . 2 2A; - v'kf 

H 7^ — r . 



3ki 



k'{2 + v') -v'kl{l-v') 



rb 9 



H ^8 
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where Ps is as defined in ( !9T]) . The zero temperature solutions and Wronskian are 



(132) 



Using ijjQ {r) = Ps/r one obtains 



P, 



P^ 



1 TlL\k^+V^kl) 



0{r- 



2iv{k^ + v^kj) 9r - 2v\l - v^)kf + kjk\-5 + llv 



9kik^L^ 



Comparing to (pOj) . one sees that 



Q 



tot 
S 



K 



,2 



7rL2[(2 + t;2)^2 



{2 + v^)k^ - v^kl{l - v^) 



v^)kl] 



(133) 



+ 



9L2 



12A; 

2v\l - v^)k\ + A;,2p(-5 + Ik;') 



(134) 



The result (1134p is more properly an expansion in rn/iL'^k), or, for v not too close to 1, 
in inverse powers of = kL'^/rn, which is the dimensionless wave- number used in p. 

To get an idea of when the next term in the expansion (11341) is significant, we have plotted 
in figure[T^,b a comparison of (I134p with numerical evaluations using the approach of [1]. The 
most immediate conclusion is that the value of K where (I134p becomes reasonably accurate 
depends significantly on v. To understand the comparison with numerics, first note that we 
are plotting not Q^* but Qe s^s defined in [1]: this quantity is proportional to minus 
the leading large k behavior in (I134p (corresponding to the pure Coulombic near-field of the 
quark) plus the non-conserved contribution from (l94l) . The constant of proportionality is the 
same as the one in (11431) . Also, because of the comparisons with away-side jet-splitting that 
motivated the study [1], we choose for an angular variable not the angle 6 = aictajik±/ki 
between the wave-vector and the direction of the quark's motion, but instead A0 = n — 6. 
Thus A(f) = TT corresponds to wave-vectors exactly along the direction of the quark's motion. 
Finally, all dimensionful quantities are rendered dimensionless using factors of ttT = th/L"^'- 
for example, K = /c/vrT. 

In figure [T}:;,d,e,f we have made corresponding comparisons of the analytic forms (11261) 
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and (11301) to numerical evaluations following [T]. As in the comparison of to Qe, there 
is a non-trivial affine relation between Qy^ and the quantity Qo appearing in the plots, and 
between and Qa- These relations can be worked out using (I143p . keeping in mind that 
the quantities plotted exclude the purely Coulombic near-field. 

5.4 The position-space near field 

The large k results from the previous sections can be Fourier transformed to position space 
without much difficulty. It is important to remember that Q^*, Qy^, and Q^* give only 
the contribution (T^n) conserved to the total stress tensor in the plasma, and that the stress 
tensor of the thermal bath {Tmn)ha.th together with a non-conserved piece (Tm„)drag that is 
responsible for the drag force need to be added in order to recover the full stress tensor 
(see equation (jTHD). From now on we will focus on the quantity {Tmn)tot = (^mn) conserved + 
(Tmn)drag, as defined in ([78]). 

Combining (I134p with the non-conserved contribution (p4l) we obtain for the energy 
density 



where we have set kf = {1 — v'^)k\ and used k"^ = k\ + k^ as in (11231) . The integrals in (11351) 
can be easily computed by differentiating appropriately modified versions of the standard 
formulas 




e{k) = e(o)(^) + e(2)(^)T2 + O(T^) 




(135) 




(136) 



The energy density in position space can thus be written 



e(t,x) 



e(°)(t,f) + e(2)(t,f)r2 + C(T^) 




(137) 
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Figure 1: Comparisons between the large K analytic approximations found in (11341) (plots 
a,b), (I130p (plots c,d), and (I126P (plots e,f) and corresponding numerical solutions of the 
linearized Einstein equations in the form studied in [1]. The unbroken colored lines come 
from numerical solutions at the values of = k/nT indicated, while the dashed lines come 
from evaluations of (11341) . suitably transformed into evaluations of Qe, Qd, and Qa as 
described in the text. The values of K used in plots c,e are the same as in a, and those in 
d,f are the same as in b. For values of K larger than the largest one shown in each plot, the 
numerical evaluations are almost indistinguishable from the K oo approximation. The 
Mach angle is indicated, but there is no reason to expect structure near this angle because 
large K is the opposite limit of where hydrodynamics is expected to apply. 
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a) e<°^ +e<^' for v=0 . 95 




Figure 2: (a) A contour plot of the near- field energy density as computed from fll37p at time 
t = 0. The red regions in the left side of the plot represent negative (Tqo) — more precisely, 
an energy deficit compared to the average energy density of the thermal bath. The deficit 
arises because the 0{T^) dominates over the 0{1) term, which leads us to question the 
accuracy of the asymptotics in these regions, (b) The (9(T^) correction e^^-* as a function of 
A(f) = TT — arctanxx/xi for several different velocities v. e*-^-* oc and a different value of 
X = x\ was chosen for each velocity in order to make the qualitative effects easy to 

discern. 

where x\ = {xi — vt)/(l — v'^Y^'^ measures the distance along the x^-axis in the co-moving 
frame, and x± = {x'^+x'^Y^'^ . At t = and for fixed v and a/ Qym^ ■> ^'^'^ ^^"^^ functions 
only of = vr — arctanx^/xi, up to overall power-law dependences on x = \/x\ + x\. In 
figure [2] we illustrate in two ways the effect of the O(T^) term on the near- field energy density. 

The main qualitative effects are Lorentz flattening of the region of large energy density, 
visible in figure [2^, and a surprising pattern of energy build-up in front of the quark with a 
corresponding deficit of energy behind the quark, visible both in figure [2^ and in the v = 0.9 
curve in figure [2)3. To understand this latter effect, recall that the subleading term e*-^^ 
encodes the leading behavior of dissipative effects at large k. For v ^ 0.67, this shift in the 
energy density is positive behind the quark (i.e. for A0 < it/2). This coincides with our 
intuition that the trailing string represents color flux spreading out behind the quark as it 
moves forward through the plasma. But for ^ 0.79, the shift in the energy density reverses 
sign as a function of 9^ A possible interpretation is that there is a "bull-dozer effect" where 

''In the intermediate range of velocities, 0.67 ^ v < 0.79, the regions where e^^^ is positive comprise a 
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the quark's near field creates a pile-up of thermal gluons and other matter just ahead of it, 
like dirt piling up in front of the blade of a bull-dozer. It would be interesting to study this 
effect further and see how far it survives beyond the leading correction term e^'^\ 

Other components of the stress tensor may also be straightforwardly computed given Q^^^, 
Qv^-i Qt^, and the non-conserved terms (l88l) . and (|95|1 . In particular, the Poynting 

vector is as follows: 



S = {Si, S2, S^) Si = 5}"^ + 5f ^T' + C(T^) 



Si{t,x) = - 
Si{k) = - 

SS) = V 

sS) = v 



tot 



(2vr) 



(0) _L 5(2)2.2 

ivkit+ik-x ^.^^^ 



V 2kl + Zk] 



v^r^^ 2Ak 



V Sniki 



V 47i{7e + 9kj 



Sikik^ 



V A7r{7k^ + 9kj 



Sikik^ 



[j.2 _^ Q^rj^A 



:i38) 



The Fourier integrals are again derivatives of the standard forms (11361) . and the position- 
space forms arc 



Si{t,x) 

S2{t,x) 

Ssit,x) 



y2 jr6 



67r2 



1 + 



TT^ XiX 2 



4 x] 



{2xi + xi)T' + 0{T'' 



X1X2 



x^ 



X 



67r2 



TT^ X 



(139) 



9ymN 



1 - 



TT^ X 



llxf + 8xl)T^ + C(T*) 



In figure [3] we show a cross-section of the approximation (I139p Poynting vector field for 
V = 0.5 and v = 0.95. For the purposes of plotting, a convenient way to regularize the 
divergence at x = is to take a cross-section at fixed but small X3. 



forward-pointing lobe and a backward-leaning, thickened cone. This is exemplified in the v = 0.7 curve in 
figure [2b- 

®We obtained the explicit position-space result (|139p for the near-field Poynting vector after seeing a 
comparable result in a summary draft of |2^. We thank A. Yarom for correspondence about this result. 
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Figure 3: The near- field Poynting vector to the approximation shown in fll39p in the plane 
ttTxs = 0.02. The arrows show the direction of the projection of S into the plane, and the 
color corresponds to the magnitude of of 5, where red means large and blue means small. 

6 Conclusions 

We have derived fully decoupled forms of the Einstein equations linearized around a general 
five-dimensional warped background ([2]) appropriate to describing holographic duals of finite 
temperature field theories on R^'^. The final equations, flMl) . fj4T|) . and (|5^ . are expressed in 
terms of diffeomorphism-invariant master fields, from which one may extract all components 
of the metric perturbations. Indeed, the entire derivation of the master equations has been 
cast in an explicitly diffeomorphism-invariant formalism. Although our results for the scalar 
master equation are incomplete in that we did not give the full expression for the source 
term except in the ArfS's-Schwarzschild background, this general source term may be readily 
obtained starting from ( l50l) . 

We have applied the master field formalism to reproduce aspects of the results of [1], 
which dealt with the stress tensor produced by an external quark moving through a thermal 
bath of A/" = 4 gauge theory. We have also extended these results by calculating the leading 
correction at large wave-number to the Coulombic near-field of the moving quark. This 
calculation was done using a method recently developed in [19]. The result for the energy 
density has an unexpected feature in position space: for v ^ 0.79, the leading correction 
indicates a pile-up of energy just in front of the quark. 

Following [T] we attempted to Fourier transform the numerical results to position space. 
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We found this difficult because even after the subtraction of the Coulombic near-field, the 
Fourier integrals have poor convergence properties: without the oscillating factor e*'^'^, they 
are quadratically divergent in the ultraviolet. With the leading correction subtracted off 
as well, this divergence is reduced to logarithmic, leading to a significantly more tractable 
numerical problem. Ideally one should calculate to one more order of accuracy in the large 
k asymptotic forms, so as to cure the divergences altogether and work with absolutely con- 
vergent Fourier integrals. Obtaining the next correction is significantly harder than the 
calculations we did in section 15. 1115. 3[ but it may be tractable. 
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A Relation to the variables used in hep-th/0607022 



In [T], the metric perturbations /i^i, were written as 



(140) 



where in axial gauge 



'-13 










Ho2 


Ho3^ 

















Hoi 





Hii 


Hl2 


Hi3 


Hq2 





Hl2 


H22 


H23 


\Ho3 





Hi3 


H23 


H33) 
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The equations of motion were decoupled into five sets of differential equations in quantities A, 
Bi, B2, C, Di, D2, El, E2, E3, and £'4, which were defined as particular linear combinations 
of the HmnS (see pQ). Setting = and k2 = k±, and decomposing the above expression 
for h^^ as in (^-(^, we get 



T,even 



y.even 



V.even 



H 



T^^n'^ k"^ 
2L^rvki 

■Di 



T,odd 



y,odd 



7-7 

^ kC 
tv — ^kikB2 



r 



H 



H 



k 



E-i 



2iL^rv ki 

' k 



Eo 



V,odd 



Jtt 



L^rv 



^ ^5 



eBi 



(142) 



H 



'0^3 



E, 



3r| 



{E3 — El] 



and of course 



y,cvcn 



y,odd 



Ji 



tr 



in axial gauge. We can use these formulas 
to relate Q^S QyS a^id Q^* to the quantities Q^*, Qd^, and Q^*, respectively, which were 
defined in [T]: 



^tot 
V5 



A;2 



k 



AkiE^v 



:2r„ ( g*/ + 



2kiL^ 



(143) 
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